Temperature effect in the conductance of hydrogen molecule 
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We present a many-body calculation for the conductance of a conducting bridge of a simple 
hydrogen molecule between Pt electrodes. The experimental results showed that the conductance 
G = dl/dV has the maximum value near the quantum unit Go ~ 2e^//i. The I — V dependence 
presents peak and dip and we consider that the electron-phonon interaction is responsible for this 
behavior. At T = there is a step in this dependence for the energy of phonons ujq which satisfies 
eV = uJo- We calculated the conductance at finite temperature and showed that dG{T)/dV oc 
l/4Tcosh2£i9^. 



I. INTRODUCTION 

The recent advances in manipulation of single 
molecules permit to measure the transport properties of 
a setup formed from an individual molecule between two 
electrodes. In contrast to the quantum dots, based on 
the semiconducting islands, the molecular devices have a 
more complicated electronic structure. 

In a remarkable experiment [l| a setup consisting from 
a single hydrogen molecule between Pt electrodes was 
measured and the conductance behavior demonstrated 
the influence of the vibrational degree of freedom in the 
transport. The effect was also observed in various or- 
ganic molecules 0, carbon nanotubes [1] fulerenes [1| 
but a theoretical description of these systems is more 
complicated because of the energy spectrum of the elec- 
trons from these molecules. The simplest model which 
is realistic for the hydrogen molecule setup is to consider 
the coupling between the vibrational mode (considered as 
phononic) of frequency lhq and electrons from the leads. 
The main effect which appears in the TV characteristic 
is the occurrence of a step at V = loq/c, which corre- 
sponds to dip versus peak in d^IJdV'^. Such a behavior 
has been studied also in d, H, Q. Many models have 
been proposed by different authors [1, H, [ll|in order to 
explain the influence of vibrational modes on the trans- 
port, but only recently Egger and Gogolin presented 
an analytical calculation explaining the current-voltage 
relation. Their calculation, which gives the correction 
to the current 51 ,is a perturbative calculation for the 
electron-phonon interaction and takes only the 
tributions similar to the approximation from Ref . | 
where the electron-phonon interaction is treated in the 
weak coupling approximation. In a higher order of per- 
turbation theory the step behavior is expected at nwo 
. The occurrence of the step feature at multiples of tJo is 
related to the strong- coupling picture using the polaron 
transformation p^ . 

In this paper we present a many- body calculation 
of the finite temperature conductance for the hydrogen 
molecule between metallic leads. The procedure is similar 
to this from Ref. [12] but we calculate the conductance 
of the setup close to the characteristic frequency &t 
finite temperature. The peak (dip)which appear at this 
frequency has a flnite width, due to thermal effect. 

The paper is structured as follows. In Sec. II we present 
the model and calculate the Green functions. In Sect. Ill 



we calculate the lowest order correction 51 to the current 
which is given by the electron -phonon processes. The 
concluding remarks are presented in Sec. IV. 



II. MODEL AND GREEN FUNCTION 

We start with the Hamiltonian which describes the 
interaction between electrons and the Holstein phonos 
used in Ref. 0] for the molecular dots and in Ref.fisj 
for study of the inelastic scattering influence on the the 
tunnelling current in the two dimensional systems. The 
simple Hamiltonian, which is analytically tractable, has 
the form : 



H = Ho+ H, 



(1) 



where 



Hq = eod^d + uob^b + ^(e^ - Ha)cl^Cka (2) 



ka 



= '^{Vad''cko, + H.c) + gQdU 



(3) 



where Q = b + b'' . 

In the Hamiltonian (1) d and are the operators for 
the single level eo, the Holstein phonons with energy loq 
are described by the operators 6^ and b, the electrons with 
the energy efeQ,(a ^ L, R) and chemical potential fj.^ are 
described by the operators Cka and c^^. The interaction 
between the electrons from the leads and the impurity 
is Va and g describes the interaction between the local- 
ized electronic level eo and the phonons. This Hamil- 
tonian has been used by Egger and Gogolin [l^] in this 
problem at T = 0. Using the equation of motion method 
we calculate the Green function Gq =<< djd'f >> as 
GUlo) = {uj ~ eo + iV)-^ where we define T = Tl + 
and Ta = 7rN{0)\Va\'^ . The Green function describing 
the system is : 



G"'(u;) = GS(cj) -I- G''(w)E''(t^)G5H 



(4) 



where the self en ergy S''(a;) is taken in the lowest order 
and has the form|14|: 
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duj'duj 



„ [/ mGo(u;', eo)Jmgo(^", c^o)]A(T, lo' , lo") 



(5) 



where A{T,uj',uj") = 1 - + /j^l^^) is the 

Fermi function and 715(0;) is the Bose function , and the 
phononic Green function Dq{u),ujq) is : 



ZuJri 



1 



1 



uj — ujQ + iS UJ + UJQ + iS 



From Eqs.(5-6)we obtain 



a,s— ±1 



Tg/f (t^O - S{fla ~ Uj)) 

{uj + SLUq) + 



(6) 



(7) 



where = {^j.r + hl)/2 - eo and Ha=L/R=±i = M±eV"/2. 

This contribution is dominant in the correction in cur- 
rent for smaU wq/F and we wiU calculate only this 
contribution at finite temperature. 



III. CURRENT AND CONDUCTANCE 

The electrical current through the dot can be calcu- 
lated from the Green function G'^{u) as : 

m = ^ / dcoihiu;) ~ fR{u)]ImG^{u) (8) 

and for the case of 5 = we get the current /o as 
e 4FrFD 

Io{V) = -^r-^[arctan(/ii/F) - arctan(/2j^/F)]. (9) 

For F — > the transparency of the junction, T = 
{h/e^)dI/dV is 



_ 4F^Ffl 1 

F 1 + {n/T f - ■ 



(10) 



Using Eqs.(4-8) we calculate the correction to the current 
given by the inelastic elect ron-phonon scattering [12] as: 



e F^Ffl 
h F 



/■Mi 

/ dujF{uj, V) 



(11) 



where F{uj) is given by : 

F2-C^2 



F{u:,V) = 



/b[^0 - s(Ma ^ ^)] 

{uj + scjo) + 



(12) 

In order to calculate the conductance G — dl/dV and its 
derivative dG/dl we will use the relation: 

d /"^^ p ff^^ dF 

Fiu;,V)^-[FmL,V) + Fmi^,V)]+J_ dcj — . 

(13) 



Using this formula we calculate the derivative d5Iinei /dV 
as : 



dV 
where : 



h F 



fF{iOo-eV)Si + e(y^-^)S: 

' (14)' 

(15) 



and 



E 



(^a - awo)^ 



i(M^ + r2)[(M^-«^o) + r2] 



212 ■ 



(16) 



Using these results we calculate the contribution of the 
inelastic scattering in ^ as : 

d^SIinei FiFij 2 



dV^ 



-9 



S1 + S2 



4Tcosh^ 



-"0 



(17) 



2T 



This is the main result of this paper, which shows that 
at finite temperature the peak (dip) in the conductance 
derivative dG/dV has a finite width at eV — ujq at low , 
but finite temperature. At T = we reobtain the 6{eV — 
ujq) behavior predicted in Ref.[12]. 

The contribution of the elastic scattering present a log- 
arithmic divergence at eV — loq which create symmetric 
dip or peak in the differential conductance. The rela- 
tive importance of the inelastic versus quasi-elastic con- 
tributions has been analyzed in at T = where was 
showed that at wq/F << 1 and large p, the inelastic chan- 
nel is dominant. This result remain valid also at finite 
temperature and in the following we will discuss the con- 
ditions for the occurrence of a dip or peak in the inelastic 
correction. First we mention that this correction is not 
any more singular, as at T = 0, but we obtain the result 
from [12] in this limit. For F^ = Tr and p. = we have 
T = 1 for Wo > 2F instead of a dip we get a peak. Thi s is 
a particular case , and for 7^ and p > ^F^ + Wq/4 we 
have a peak. We also have at p, = ±F the transparency 
T = 1/2 which is in fact the point of the peak-dip tran- 
sition. 

At T = the only way to obtain a finite life-time of 
phonons , which may give a smearing of the step and/or 
peak feature is to include the electronic polarization in 
the Green function of the phonons, as was suggested in 
However, such a calculation implies the higher-order 
perturbations in g and it is difficult to be performed an- 
alytically. In this paper the discussion concerning the 
occurrence of the step or/and peak is identically to that 
of the authors of [T^, but the coupling to the thermal 
phonons generate in a natural way a smearing of the step 
or/and peak structure. 
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IV. CONCLUDING REMARKS 

We analyzed , using the many body method, the trans- 
port in the simplest molecular dot consisting from a hy- 
drogen molecule between the Pt leads. The analytical 
calculations of the effect given by the electron-phonon 
interaction on current at finite temperature have been 
performed. The obtained results can be regarded as com- 
plementary to the T ~ similar calculations presented 
in Rcf. 12]. However, at zero temperature the concept 
of phonon is not defined and in the transport we use 
the concept of inelastic scattering. On the other hand the 
experimental results showed the existence of behavior 
in the dG/dl which cannot be described by the simple 
6{eV — ujq) behavior, which is specific for T — 0. We 
showed that at finite temperature a dip or a peak in this 



quantity is described by dG/dl = C{fl,T,u!o)l/{4T 
cosh?' ^^~r^° ).This behavior is given by the inelastic scat- 
tering between electrons and phonons, the elastic contri- 
bution giving a non relevant contribution to the trans- 
port in this system, which can be considered as a molec- 
ular quantum dot. Our results at T = are identically 
with the results from [T3|.and the signature the constant 
C(/2,r,wo) gives the same conditions. Our calculations 
completed the microscopic model presented in Ref.jTl], 
but we consider that the transport in more complicated 
molecules is difficult to be treated analytically.However, 
the model can be a starting point for the study of trans- 
port in the more complex molecular systems. 

We thank Alexander Gogolin for useful correspondence 
on the subject. 
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